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We study the relaxation of the bi-dimensional kinetically constrained spiral model. We show 
that due to the reversibility of the dynamic rules any unblocked state fully decorrelates in finite 
, times irrespectively of the system being in the unjammed or the jammed phase. In consequence, the 

O^l ■ evolution of any unblocked configuration occurs in a different sector of phase space from the one that 

includes the equilibrium blocked equilibrium configurations at criticality and in the jammed phase. 
We argue that such out of equilibrium dynamics share many points in common with coarsening in 
the one-dimensional Ising model and we identify the coarsening structures that are, basically, lines 
• of vacancies. We provide evidence for this claim by analyzing the behaviour of several observables 

including the density of particles and vacancies, the spatial correlation function, the time-dependent 
persistence and the linear response. 
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I. INTRODUCTION 



Kinetically constrained models are toy models for the glassy phenomenon [l[ (see for reviews). These models 

display no thermodynamic singularity: their equilibrium measure is simply the Boltzmann factor of independent 
variables and correlations only reflect the hard core constraint. Still, they capture many features of real glass- forming 
systems; among them one can mention stretched exponential relaxation [f|, super- Arrhenius equilibration time [||, 
dynamical heterogeneity Q and self-diffusion/ viscosity decoupling They also display glassy out of equilibrium 
aspects, such as physical aging Q and effective temperatures [lO] (for a review of these aspects see that are similar 
t~^. ' to the ones first obtained in random p-spin systems [IJ EH > as weu as heterogeneous aging relaxation [l3| . Bootstrap 
, percolation arguments allowed C. Toninelli et al. to show that most finite dimensional kinetically constrained models 
studied so far do not even have a dynamic transition at a particle density that is less than unity in the thermodynamic 
limit [14]. On Bethe lattices instead a dynamical transition similar to the one predicted by the mode coupling theory 
f-*. ■ might occur (l5j . 

A finite dimensional kinetically constrained model with an ideal glass-jamming dynamic transition at a particle 
density that is different from one, the two-dimensional spiral model, was introduced in a series of papers by C. 
Toninelli et al. (l6l . UtI [l8l fl9j | . The dynamic transition is defined as the critical density at which an equilibrium 
configuration can no longer be emptied by applying the dynamic rule. In consequence, at the transition an infinite 
cluster of particles (or spins in an equivalent representation) that are blocked by the dynamic constraints exists. The 
model is defined on a square lattice and the transition occurs at a value of the control parameter p, the probability 
to update an unblocked site, that coincides with the critical threshold of directed site percolation in d — 2, that is 
to say, p c ~ 0.705. The density of the frozen cluster is discontinuous at p c which means that the frozen structure 
is compact rather than fractal at criticality. Thus, at and above p c a finite fraction of particles are blocked. It was 
shown in these papers that approaching the critical p c from the unjammed phase the decorrelation of an equilibrium 
configuration decays in two steps similarly to what is common in glass- forming liquids. The time-scale for relaxation 
- say the a relaxation time - diverges with a Vogel-Fulcher law (faster than a power law). A crossover length below 
which finite size effects are important also diverges faster than a power law when p c is approached from below. 

In this paper we study the dynamics of the 2d spiral model numerically. For concreteness, we first show results for 
the equilibrium dynamics, meaning the relaxation of an equilibrium configuration, below but close to the critical p c - 
the 'super-cooled liquid' - and above p c - the 'equilibrium glass'. Next, we focus on the out of equilibrium kinetics. 
Specifically, we consider annealing and quenching processes. In the former the system is prepared in an equilibrium 
initial condition at po > p c and subsequently brought into the glassy regime with p > po by using a finite rate of 
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FIG. 1: (Color online.) The neighbouring sites determining the frozen or free to move character of the center site 

change of the parameter p. In the latter we study the evolution after a sudden quench, in which an equilibrium state 
at po = (completely empty configuration) is evolved from time t = Q onwards with the dynamic rule specified by a 
different value of p > p$. This procedure is similar to the temperature quench of a liquid. For p < p c , after a non 
equilibrium transient the system attains the unblocked equilibrium state. The relaxation time for attaining such a 
state diverges as p — > p c . In the critical case with p — p c , the system approaches the blocked equilibrium state by 
means of an aging coarsening dynamics with some points in common with the ones observed in critical quenches of 
ferromagnetic models. Freezing is not observed because the dynamic density p(t) is always smaller than the critical 
one p c — Pc, at any finite time. A different situation occurs for filling at p > p c . Also in this case the system keeps 
evolving up to the longest observed times and a blocked state is never observed despite the fact that p(t) exceeds p c 
for long enough times. This can be understood by noticing that the dynamics are fully reversible: for any p, each 
configuration reached dynamically can always evolve back to the initial state with p < p c by the time reversed process, 
although with a very low probability. Therefore, states with blocked regions cannot be dynamically connected to the 
initial unblocked states and they cannot be reached by any reversible process. Frozen regions being dynamically 
forbidden, the system approaches a high density (p > p c ) ensemble which does not contain any of the blocked states 
the measure of which is relevant in equilibrium. This feature recalls what happens in coarsening systems where, again 
for dynamical reasons, the kinetics occurs on the ensemble of configurations with vanishing magnetization the weight 
of which is negligible in equilibrium [20] . Actually, it will be shown that the dynamics of the spiral model at p > p c 
strongly resembles coarsening in ferromagnets. 

II. MODEL 

A binary variable = 0, 1 is defined on the sites of an C x C square lattice in d — 2. = 1 corresponds 
to a particle, riij — to a vacancy. For a given position let us define the following couple of neighbouring sites 

(see Fig. [TJ 

• + + l,j + 1), north east (NE) couple 

• (i + + 1), south east (ES) couple 

• — — 1), south west (SW) couple 

• (i — — l,j + 1), north east (WN) couple 

The possibility to update the variable depends on the configuration of such couples. Specifically, if the sites 
belonging to at least two consecutive couples (namely, NE+SE or SE+SW or SW+NW or NW+NE) are empty 
(otherwise stated occupied by vacancies), then site can be updated (either emptied or filled). Otherwise it is 
blocked. We consider each site coupled to a particle reservoir in such a way that particles can enter or leave the 
sample from its full volume, in contrast to what was used in Q, for instance, where particles were allowed to access 
or leave the (3d in this case) system only through the borders. We use periodic boundary conditions. 
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FIG. 2: (Colour online.) Equilibrium configurations at (a) p = 0.5 < p c , (b) p c ~ 0.705, and (c) p = 0.99 > p c . The colour 
code is such that: a filled blocked site is black; an empty blocked site is white; a free particle (one that can be removed) is red 
and a free vacancy (one that can be filled) is green. We show a small fraction of the sample to make the snapshot easier to 
visualize. 

Defining the variable My- such that My = if site is frozen, and My = 1 otherwise the updating rule can be 
compactly expressed in terms of the transition rates Wp(ny|ny) to go from n to n' on site as 

W p (0\l) = My p W p (l\0) = My (1 - p). (1) 

The equilibrium state has been studied in detail by C. Toninelli et al. [H, The model can be regarded as a 
two-level system described by the Hamiltonian H [s] — — y\ . ny at the inverse temperature 

/? = ln[p/(l-p)] . (2) 

When p — ► 1/2 the inverse temperature vanishes /3 — -> whereas for p — > 1 it diverges /3 — > oo. The Bernoulli measure 
implies p = (ray) = p for the equilibrium density of particles. This means that, in equilibrium at zero temperature 
(p — > 1), the lattice is full while at infinite temperature (p — ► 1/2) it is half-filled. The existence of a blocked cluster 
at p > p c < 1 was shown by studying the critical p above which an equilibrium configuration cannot be emptied. 

Interestingly enough the microscopic dynamic rules are fully reversible. This implies that an initial configuration 
with no blocked structure can never reach, via the dynamics, a state with a blocked sub-ensemble of sites within. This 
fact holds for a system coupled to the reservoir in any possible way (volume, border, or other). We shall elaborate on 
the consequences of this fact later in the paper. 

III. OBSERVABLES 

The dynamics of kinetically constraint models are usually studied in terms of the probability that a spin never flips 
between time t w and a later time t, also called the persistence function <f>(t, t w ). Operatively, we define this quantity 
as 

<j>{t,t w ) = —(N noup (t, t w )), (3) 

where N — C x C is the number of lattice points, N noup {t,t w ) is the number of sites that have not been updated 
between t w and t in a given realization of the dynamics and (...) denotes an average over different histories and 
initial conditions. In equilibrium conditions one has (j)(t,t w ) — <fi(t — t w ). Moreover, at and above p c one has 
= limt^oo (f>(t, t w ) > due to the frozen backbone, while below p c the asymptotic value vanishes </>oo = 0. 
Similar properties are exhibited by the autocorrelation function 

C(t,t w ) = (nij(t)riij(t w )) - (nij(t)){n,ij(t w )}. (4) 

The impulsive autoresponse function is defined as 



eij=0 



(5) 
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where e,j is a perturbation changing the transition rates as W p (nij , n'^) — > W p + eij (riij , nj,-). Notice that we dropped 
the site dependence in C and i? due to space homogeneity. The integrated response function, or susceptibility, is 

X (t,t w ) = f dt'R(t,t') . (6) 



In equilibrium conditions the response function is related to the autocorrelation function by the fluctuation-dissipation 
theorem (FDT) 

T X (t,t w )=C(t,t)-C(t,t w ) (7) 

where T — with f3 defined in Eq. ^fy in the spiral model. Out of equilibrium the response function can still be 
related to correlation functions of the unperturbed system (with e = 0) according to the generalization of the FDT 
derived in [22] . This will allow us to determine the response function numerically without applying any perturbation. 
It is useful to introduce the normalized autocorrelation function 

C(t,t w )=C{t,t w )/C(t,t) (8) 

and susceptibility 

X(t,t w )= x (t,t w )/C(t,t). (9) 

Finally, let us define 

l(t) = [p eg -p(t)}- 1 . (10) 

In the non-equilibrium cases considered in Sec. IV Bi when an initially empty state is progressively filled with the 
transition rates W p in Eq. ([I]), p(t) always approaches p eq = p. In this case l(t) is a monotonically growing function 
of time, and it will be convenient to re-parametrize t in terms of I. 

The spatial evolution of the structure is usually monitored via the equal time correlation function (equivalently, in 
Fourier space, the structure factor): 

G(r,i) = ( nij (t)n iff (t)) ~ {n t] (t)) 2 (11) 
where r is the distance between sites and (i',f). A typical length can be extracted from G by using 



L(t) 



f dr r a G(r, t) 



-i 1/a 



J dr G(r, t) 



(12) 



When scaling holds, namely G(r, t)/G(0, t) = g(r/L(t)), changing a leads to the same determination of L(t) apart 
from a multiplicative constant. 

All our simulations are done with square systems with linear size C — 200 — 2000, depending on the situation. We 
typically average over 10 2 — 10 4 realizations of the dynamics. 



IV. EQUILIBRIUM RELAXATION 

Figure [2] displays three panels with equilibrium configurations at p < p c (a), p = p c (b) and p > p c (c). We shall 
later confront these images to the ones reached dynamically from an empty state, see Fig. [8] 

In Fig.[3{a) we study the relaxation of an equilibrium configuration for several values of p by means of the persistence 
function. For values of p < p c there is no frozen structure in the system and the persistence function, (f>, decays to 
zero. Close to the critical p Cl although a blocked spanning structure does not exist, long living quasi-frozen structures 
do. The decay of the persistence or the normalized correlation close (but below) criticality is reminiscent of the 
hallmark of the super-cooled liquid relaxation that is associated to the formation of long-lived cages surrounding each 
particle (23j. The persistence of these structures diverges as p c is approached. As a consequence, <p(t,t w ) develops a 
plateau that becomes longer and longer when p gets closer to p c as in a super-cooled liquid. This numerical result 

confirms the analytic ones in [H [13, IS EH. 

Above p c the initial configurations have a finite fraction of frozen particles so <f> decays to a constant, 0oo, that 
depends on p. Figure [3fb) shows 4>oo(p). As expected (poo is different from zero at p c ; we found (/>ooipc) — 0.9, a 
relatively large value. A fit to the data close to the critical point yields (f>oo{p) — <t>oo(Pc) + c(p — p c ) a with a ~ 0.9. 
increases with p and reaches 1 at p — > 1. 

The normalized autocorrelation function C(t,t w ) behaves in a qualitatively similar way. Moreover, consistently 
with equilibrium dynamics C(t,t w ) is related to x(t,tw) by the FDT, namely T X (t,t w ) = 1 — C(t,t w ), see Eq. ([7]). 
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FIG. 3: (Color online.) Upper panel: Decay of the persistence function for an initial condition in equilibrium at different values 
of p. 4>oo = linit^oo 4>(t) = below p c and linit^oo if)(t) > above p c . Lower panel: Dependence of 0oo on p. In the inset 
4>ca(p) — 0oo (p c ) is plotted against p — p c in a double logarithmic scale, for p > p c . The red dashed line is the best fit to the 
power-law behavior 4>oa(p) — <f>oo(Pc) oc (p — p c ) a , yielding a — 0.901. 



V. OUT OF EQUILIBRIUM RELAXATION 



In this Section we discuss the out of equilibrium dynamics of a system initially prepared in some equilibrium 
configuration at pq and then evolved from time t = onward with the transition rates |T]) and p > po- As discussed 
in Sec. [II] this procedure is the analogue of a cooling experiment in thermal systems. 

We shall consider two types of 'cooling' procedures: the former is a slow process where a system at equilibrium at 
p = Po > Pc, a configuration with a blocked structure, is evolved with the transition rates (|TJ) and a time-dependent 
V = Po + ft, r being the equivalent of a cooling rate. The latter is an instantaneous quench where an equilibrium 
system at po = 0, a completely empty state, is evolved using p < p Cl p = p c and p > p c as target values of the control 
parameter for the subsequent dynamics. 
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FIG. 4: (Color online.) The density of particles as a function of p for linear cooling rate procedures pit) = po + rt starting 
from equilibrium at po > p c - Continuous curves: po — p c — 0.705 and different cooling rates r given in the key. The annealed 
densities at fixed p increase with decreasing r but they remain well below the equilibrium values, p eq = p, shown with an orange 
straight line in the figure. Dashed line: the initial condition is equilibrium at po = 0.71 > p c and the cooling rate is r — 10~ r . 
In the inset the continuous and the dashed line represent the curves obtained with r — 10 -7 starting from po = p c (continuous 
line) or p = 0.7051 > p c . 



We recall that the dynamic rule is reversible and configurations with or without a blocked sub-ensemble of sites are 
mutually inaccessible even for finite systems, a fact that draws an important difference with usual stochastic models 
without dynamical constraints. 



A. Slow cooling from equilibrium at p c : threshold level 

In Fig. 2] we show the density of particles found after cooling with a finite rate. Continuous lines with different 
colour represent processes starting from an equilibrium configuration at p c and using different cooling rates r. Dashed 
curves, instead, refer to cases where cooling starts from equilibrium above p c . Concerning cooling from equilibrium 
at p c , for p's that are above but not too far from p c the curves for r < 10~ 7 , say, superpose. For larger values of 
p there is a residual r dependence that presumably disappears for even smaller r's. We assume then that the dark 
green curve is quite close to the one for the limit of vanishing cooling rate pth(p) = lim r ^o Prip)- Pthip) is definitely 
below the equilibrium density p eq for all p > p c and hits a limit pthiX) ~ 0.723 < 1 when p — » 1. 

This analysis demonstrates that the blocked structure at p c involves vacancies that cannot be filled when increasing p 
with the cooling procedure and the dynamic density remains well below the equilibrium one (shown with a continuous 
straight orange line in Fig. [J). A similar procedure in which T is changed infinitely slowly, carried out in the p-spin 
model with p > 3 defines the threshold free-energy level [llT ] that remains higher than the equilibrium free-energy in 
the full low temperature phase. 

Cooling with a finite rate an initial condition equilibrated at a value po(> p c ) yields other levels of blocked states 
(low lying metastable TAP states in the p > 3-spin analogy) . In Fig Q] we show with a dashed line one of such curves 
with po = 0.71. A more careful analysis shown in the inset suggests that the curves for different po do not cross in 
the full p > p c region, again in accordance with the p-spin scenario in which there is neither level crossing nor level 
merging when T is lowered infinitesimally slowly. 
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FIG. 5: (Color online.) The difference between p(t) and its equilibrium value p eq against time for different values of p. The 
relaxation is faster the smaller the value of p. 

B. Infinitely rapid quench from equilibrium at po = 0. 

In this Section we consider the dynamical process whereby an empty configuration is filled with the transition 
rate ([T]) , for different values of p. We have checked that starting from an equilibrium configuration corresponding to 
< po < p c (i.e. not completely empty), does not change the evolution significativcly. In the case of quenches to 
p > p c the effect of po > is to delay the asymptotic regime where, as we shall discuss below, scaling holds. This 
is due to the more blocked nature of the initial state. For this reason, the choiche po = is the most convenient to 
study the asymptotic properties, and we shall always make this choice in the following. 

In Fig. [5] we plot the difference between p eq and the dynamic density pit). This quantity is related to l(t) through 
Eq. (|10p . For p < p c the density decays exponentially to the equilibrium value over a typical time that increases as p 
approaches p c . At p = p c , there is a power law decay p(i) — p eq oc t~ v ( pc ^ with v(p c ) — 0.16. For p > p c , the same 
behavior, p(t) — p eq oc t~ v ( p > is observed, although the exponent vip) becomes smaller as p increases and seems to go 
to zero in the limit p — ► 1. Notice that this implies that lit) grows as a power law in the whole phase p > p c , and 
logarithmically in the extreme limit p — ► 1. The relaxation seems to be smooth in the sense that p(t) decays as a 
power law within the whole explored time-window. 

If one extrapolates the power law behaviour to infinitely long times the fact that p eq — pit) — > would mean 
that, surprisingly enough, the relaxation avoids the threshold level - defined by the procedure described in Sect. IV A"l 
Part of the explanation of this fact is that the dynamics are reversible and after a quench the system never reaches 
configurations with a completely blocked sub-ensemble. This implies that the system cannot go through threshold 
states to lower lying metastable ones since all the former do have blocked sub-ensembles. This is clearly different 
from what is the p-spin scenario in which a typical initial condition as the one corresponding to infinite temperature 
first relaxes to the threshold level [TJ HH H3] and then penetrates below it via thermal activation for finite size 
systems. It is also different from the behaviour of kinetically constrained spin models on the Bethe lattice [lj| that 
seem to conform to the p-spin scenario. The fact that the out of equilibrium relaxation of the spiral model reaches 
pit) > pthip) implies that the dynamics should go to the long-time configurations 'turning around' blocked ones in 
phase space. 

In the following Sections we shall discuss the three cases of quenches to p < p c , p — p c and p > p c in more detail 
separately. 

1. Quench from po = to p < p c : equilibrium dynamics. 

We have already observed that for p < p c but sufficiently close to p c long lasting quasi-blocked structures exist 
which are responsible for the development of a plateau in the persistence function. In the case considered here, since 
the initial configuration is empty, those structures must be built. Hence the plateau observed in equilibrium (Fig. [3]) 
can be observed only in the limit t w — > oo, when equilibrium is attained and the long lived quasi-frozen structures are 
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FIG. 6: (Color online.) The persistence function cj}(t,t w ) is plotted against t — t w in the case of a quench from an empty 
configuration to p = 0.65, for different values of tw The equilibrium curve (long dashed) is approached in the large-t™ limit. 
The system is never blocked and <j> decreases to zero in all cases. 



fully developed. For finite t w , on the other hand, those structures are not completed and the plateau has a shorter 
extent, as shown in Fig. [6l C(t,t w ) behaves similarly. 



2. Quench from p — to p = p c : Critical quench. 



In this case a configuration with a blocked structure is approached in an infinite time but never reached. Since, 
as shown in Fig. [5l p(t) approaches the equilibrium value asymptotically as a power law, we shall re-parametrize t 
in terms of l(t) in the following. We do this because, although t and I are trivially related in the asymptotic time 
domain, at early times when corrections to the power law behavior of I are present, this quantity turns out to be 
better suited to enlighten the scaling properties of the dynamics. Since (f)(t, t w ) and C(t, t w ) display the same scaling 
properties we discuss in the following the behavior of the latter quantity. The data of Fig. [7] show that C(t, t w ) obeys 
the scaling 



C(t,t w ) = f 



l(t w ) 



(13) 



The collapse of the curves for different t w is very good for large values of x — l(t)/l{t w ) ~ 1.75, while it is poorer for 
smaller values of x perhaps due to preasymptotic effects. The same scaling is found for the susceptibility, as testified 
by the fact that the parametric plot of xit^w) against C(t,t w ) (upper inset of Fig. [7]) does not depend on t w . The 
parametric plot has the linear equilibrium behavior T c x(t,t w ) = 1 — C(t, t w ), as in ferromagnetic models quenched 
to T c . However, while in the latter case the quantity X(t,t w ) = T R(t,t w ) / {dC / dt w ) is a non trivial function oit/t w 
with the limiting value lim^^oo lim^oo X (t, t w ) = Xoo < 1 p7| . here X(t, t w ) slowly approaches X(t, t w ) = 1 in the 
limit t w — » oo, as shown in the lower inset of Fig. [7J 



3. Quench from po = to p > p c : Coarsening. 

We now discuss the dynamics following an infinitely rapid quench from equilibrium at po — to p > p c . We argue 
that the relaxation occurs through a coarsening process which resembles some features of ferromagnetic systems. The 
kinetics simplifies in the limit p — > 1 because Wi(l|0)/W r i(0|l) — » and the annihilation of particles can be neglected 
with respect to creation events if the latter are available. We shall consider this case first, since a cleaner description 
of the kinetics and of its basic mechanisms is possible. The case with a non-vanishing 1 — p finite, where the finite 
annihilation probability adds new mechanisms to the kinetics, will be considered further below. 

Quenches to p — > 1 
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FIG. 7: (Color online.) Dynamics after a quench from po < p c to p c . C(t,t m ) is plotted against l(t)/l(t m ). The upper 
inset shows the fluctuation-dissipation plot of x(t,t m ) against C(t,t w ). The lower inset shows the fluctuation-dissipation ratio 
X(t,t w ) as a function of l(t w )/l(t). 
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FIG. 8: (Color online.) Left and central panel: A zoom on the configuration of the system at two times (t = 10 5 and t — 10 r ) 
after a quench to p — 0.99. Black and white sites are frozen particles and vacancies, respectively. Red and green sites are 
particles and holes that can be updated. Right panel: Configuration of the system at time t = 10 7 after a quench to p = 0.85. 
The structure is fuzzier than in the p = 0.99 case. 

In this section we study the p — ► 1 limit by letting p — 0.99 in the numerical simulations. In the limit p — > 1, p(t) still 
decays to zero, see Fig. [5] Hence, for large times there is a very small fraction of vacancies. Despite this fact, vacancies 
cannot be isolated, otherwise the dynamical constraints would freeze the system, which is not observed at any time 
in our simulations. These two requirements can be both fulfilled if vacancies segregate into quasi-unidimensional 
domains (see Fig. El) • A representation of a quasi-unidimensional domain is given in Fig. [5] Such a domain is basically 
a set of vacancies each of which is surrounded by n = 1, 2 neighbouring vacancies. The word quasi refers to the fact 
that these strings of vacancies can be seldom decorated by some bi-dimensional feature such as a traveling vacancy 
(TV), and T-junctions, as shown in Fig. [5] Since these objects are basically one-dimensional the system is allowed to 
build very long strings with a vanishing vacancy density in the thermodynamic limit. 

In the limit p — ► 1 the dynamics proceed as follows: starting from a generic configuration at large times the only 
possible move is almost always the removal of a non frozen particle. This process takes a typical time tq ~ (I — 
which diverges as p — > 1. There is basically no evolution on time-scales shorter than To which must be regarded, 
therefore, as the shortest time in the system. Whenever a particle is removed it is always possible to add a particle (in 
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the site from which it was previously removed or in the neighbourhood) and this event almost always occurs for p — > 1. 
On the other hand, processes such as the consecutive removal of two (or more) particles, which have probability of 
order (1 — p) 2 (or higher powers of 1 — p) can be discarded. Therefore we can focus on the leading events to lowest 
order in 1 — p and the dynamics are provided by the following mechanisms (see Fig. [§]) 

• Diffusion and annihilation of traveling vacancies: with probability 1 — p one of the particles surrounding a 
traveling vacancy can be removed. When this occurs a pair of adjacent TVs is formed and, since p ~ 1, one of 
the two is immediately filled with a particle. This process, which taken alone does not change p, makes TVs 
diffuse along flat parts of the string. If two TVs come together one of them is immediately filled with a particle 
and p decreases. 

• Deformation around kinks or T-junctions: a particle can be removed with probability 1 — p from the internal 
part of a kink or a T-junction. Hence the neighbouring site can be filled by a particle. This process produces 
the deformation of the string, since kinks and T-junctions can diffuse and/or generate new kinks and TVs. In 
addition, if two kinks or T-junctions meet, a new particle can be added lowering the density. Notice that these 
mechanisms, besides lowering the density, may give rise to the coalescence of strings, increasing the size of the 
quasi-unidimensional aggregates of vacancies, a phenomenon which is clearly observed in Fig. [5] 






FIG. 9: (Color online.) Schematic evolution of the system. Black sites are particles, (a red circle indicates the ones which are 
not frozen) white sites are vacancies (a green cross indicates those which are not frozen) . Left panel: the mechanism producing 
the diffusion and annihilation of two TVs. Central panel: the mechanism producing the diffusion and coalescence of two kinks. 
Right panel: two T-junctions can diffuse similarly to kinks (see central panel) and merge. 

As discussed above the deformation and coalescence between strings lead to a perpetual growth of the typical 
string length, and consequently of the average distance between strings L(t), as can be observed in Fig. [8l This 
suggests that dynamical scaling may set in for large times, due to the prevalence of L(t) over any other length, 
similarly to what happens in coarsening systems. According to the scaling hypothesis, configurations of the system 
look statistically similar if distances are measured in units of L(t). For the equal time correlation function this 
implies G(r,t)/G(0,t) — g(r/L(i)). Similarly, the autocorrelation function must depend on t,t w only through the 
ratio x = L(t)/ L(t w ), namely 



C(t,t w ) = f(x) 



with 



Lit) 
L{t w ) 



(14) 



In the following we check this hypothesis. 

Due to the one-dimensional geometry of the strings, we expect the typical distance between them to be proportional 
to the inverse of the vacancy density p v (t). Since in this limit p v (t) = 1 — p(t) ~ ^ _1 (i) we claim 



L(t) oc l(t). 



(15) 



This is demonstrated in Fig. ITT7I where we plot L(t) computed through Eq. (fT!?)) for different values of a against l(t). 
In the presence of scaling, as already observed in Sec. IIII1 determinations of L(t) using different values of a give 
the same result, apart from a multiplicative constant. Notice however that, in the presence of noise, the quality of 
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FIG. 10: (Color online.) The typical length defined in Eq. (|12|l . for different values of a given in the key, against l(t). 

the determination of L(t) may strongly depend on the value of a. Actually, different values of a weight differently 
regions with different r and, since data at small r are typically less noisy, smaller values of a lead to cleaner results. 
This is clearly observed in Fig. 1101 On the other hand, if preasymptotic corrections to scaling are present, deviations 
from scaling can be enhanced or suppressed by changing a, because those corrections may be more effective in the 
large or small r regions, depending on the system. Actually, the data for a — 1 and a — 0.5 show no appreciable 
deviations from the asymptotic law L(t) oc l(t) in the whole time range; with a — 0.25 one observes small deviations 
up to l(t) ~ 4 — 4.5 which can be interpreted as due to preasymptotic corrections to scaling. In any case, data for 
l(t) > 4.5 clearly show that L{t) oc l{t) for all a. This implies that scaling is obeyed and that, in this case, lit) has 
the meaning of a length. 

Since G(r, t) is too noisy to obtain further evidence for scaling we consider in the following two-time quantities. In 
the inset of Fig. [11] we plot C(t,t w ) against t — t w for different values of t w . One clearly observes an aging behavior. 
The behaviour is very similar to the one in the Id Ising chain. The apparent plateau slowly moves upwards and tends 
to 1. In the main figure C(t, t w ) is plotted against l(t)/l(t w ). Although the data collapse does not provide a definitive 
evidence for scaling, there is a clear improvement upon the collapse as t w increases, suggesting that scaling may be 
asymptotically obeyed for very large t w . Interestingly all the curves intersect at the point of coordinate x ~ 1.08 with 
curves for longer (shorter) t w lying below (above) curves for shorter (longer) t w before (after) the crossing point. This 
is somewhat similar to what was found in Monte Carlo simulations of the Sherrington-Kirkpatrick (SK) spin-glass ^29] , 
that is known to satisfy dynamic ultrametricity asymptotically [301 ] . However, in the spiral model the distance between 
consecutive curves reduces as t w increases (contrary to what is seen in the SK case) thus suggesting that this peculiar 
fact is just a pre-asymptotic effect. 

Next we consider the response function. In Fig. [T2] the fluctuation dissipation plot Tx(t,t w ) against C(t,t w ) is 
shown. The curves collapse quite well on a mastercurve X P ~i(C). This implies that the response function obeys the 
same scaling as C(t,t w ), namely x(t,t w ) — h(x) with x = L(i)/L(t w ). In the short time regime, for values of the 
x axis that are close to one, one has x(C%~i ~ 1 — C. As C decreases the curve bends. Although, naively, one 
could think that the parametric plot is one with two straight lines (one with the temperature of the bath and another 
one with a different value with a sharp crossover in between) the interpretation of the dynamics in terms of one 
dimensional coarsening suggests instead that the asymptotic construction is a non-trivial curve analogue to the one 
found in the Id Ising model. Actually, as shown in Fig. 1121 Xp~i(C) is quite similar to the parametric plot xt~o(C) 
(known analytically [251 ]) of \ vs C i n the Id Ising model quenched to T ~ (no normalization is needed in this case 
since C(t,t) = 1). Since xt~o(C) is known not to be universal (depending, for instance on the observable used to 
compute the response [281 ]) we do not expect the curves for the two models to be quantitatively equal, but only to 
behave in a qualitatively similar way. On the other hand, since the limiting slope = dxT~o{C) / 'dC\c=o = 1/2 
is expected to be universal 0, [27| , it would be interesting to compare this value with the limiting slope of the spiral 
model. 



Quenches to p c < p < 1 
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FIG. 11: (Color online.) The system is prepared in equilibrium at pa = and then instantaneously quenched to p = 0.99 at 
t = 0. C(t,t w ) is plotted against l(t)/l(t m ) for different values of t w given in the key. In the inset, the continuous lines are 
the same curves C(t,t w ) of the main figure but plotted against t — t w . The dotted lines represent C(t,t w ) against t — t w (for 
different values of t w given by the same colour code of the continuous lines) after a quench to p = 0.85. 

In this Section we consider quenches from an empty state to p c < p < 1. With respect to the case p — > 1 considered 
above, there is now a non-vanishing annihilation probability, due to a finite 1 — p. Hence, after a time of order 
r n ~ (1 — p)~ n , processes where n vacancies are created start to be observed. Then, for t < T2 the dynamics proceeds 
similarly to the case of a quench to p — > 1; later on a bubble of empty sites can be created around a quasi-dimensional 
structure, as can be seen in the right panel of Fig. [8j These processes produce a local rejuvenation of the system, 
since in the bubbles the density decreases (although globally it is still growing as shown in Fig. As an effect, for 
a given t w , the autocorrelation function decays faster than in the case with p — ► 1 , as can be seen in the inset of Fig. 

HU 

The process of bubble formation can be compared to the fast quasi-equilibrium processes occurring in ferromagnetic 
systems when quenched to a finite T > 0. In scalar models (e.g. the Ising model) these processes are basically 
the reversal of spins in the bulk of ordered domains due to thermal fluctuations. In systems with a finite critical 
temperature T c (say, the Ising model in d > 1) quenched below T c , in the large-i^ limit the timescale of fast processes 
is widely separated from the timescale over which irreversible or aging events occur. Quasi-equilibrium processes are 
responsible for the fast initial decay of C(t,t w ) to a plateau value C(t,t w ) = m 2 , namely the squared equilibrium 
magnetization, which occurs for t — t w < T eq , where r eq is the equilibrium relaxation time. Being equilibrium 
fluctuations in nature, these processes obey the FDT, and Eq. ([7]) is observed in a restricted time domain when 
C(t,t w ) > m 2 . Despite this fact, the system is globally aging, as shown by the further decay of C(t,t w ) from the 
plateau value on much larger timescales t — t w > t w 3> r eq when aging processes become relevant. Notice that, by 
raising T towards T c the region of the parametric plot where the FDT is observed extends to lower and lower values 
of C until at T c Eq. ((TJ) is obeyed in the whole range of C variation. For system with T c = (i.e. the Id Ising model) 
quenched to a low T > 0, instead, the appearance of thermal fluctuations in the bulk of domains at a certain time t eq 
is accompanied by the global equilibration of the system, so that aging is interrupted. Eq. |[7J) is then globally obeyed 
from t w ~ t eq onwards. 

In the spiral model with p < 1 bubble formation seems to play a role broadly similar to thermal fluctuations in 
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FIG. 12: (Color online.) Left panel: Tx(t, t w ) against C(t, t w ) for different values of t w , in the case of an infinitely rapid quench 
from equilibrium at po = to p = 0.99. The bold purple line is the equilibrium behavior T\(t, t w ) = 1 — C(t,t w ). The dashed 
blue line is the parametric plot of x vs C m the Id Ising model. Right panel: Tx(t,t w ) against C(t,t w ) for t w = 10 3 in the 
case of an infinitely rapid quench from equilibrium at po — to different values of p > p c given in the key. 

ferromagnets, although combining features of systems with T c — to those of the cases with T c > 0, and with some 
notable differences. Actually, when bubbles begin to be formed C(t,t w ) starts to decay faster than in the case with 
p — > 1, but we could not detect the existence of a two step relaxation (with a plateau in between) in our data for 
C(t, t w ) (see Fig. [II]). This suggests that an effective time-scale separation between fast and slow processes as observed 
in ferromagnetic systems is not present in the spiral model, at least in the time sector studied in our simulations. In 
Fig. [12]- right the parametric plot of x against C is shown for quenches to different values of p > p c . Only the value 
t w = 10 3 is shown, we found that for each choice of p the curves with values of t w up to t w = 10 5 are indistinguishable 
within numerical errors from the corresponding ones with t w — 10 3 . Although we cannot exclude that a tiny in- 
dependence could be observed by considering much larger values of t w , this strongly suggests that the curves with 
different t w collapse on a p-dependent mastercurve Xp(C'). Notice that by lowering p towards p c , X P {C) rises towards 
the FDT line ([7]) and the region of large C where the FDT is approximatively obeyed is enlarged. This is partly 
similar to what observed in ferromagnets with T c > quenched below T c when T is raised towards T c . However, 
at variance with those systems, the shape of the plot is non trivial (i.e. all the curves bend continuously) and there 
is no apparent dependence on t w , although from the inspection of C one concludes that the system is still globally 
aging. These features resemble what observed in ferromagnets with T c = 0; in the spiral model however this behavior 
is observed in the whole glassy phase p > p c . 

VI. CONCLUSIONS 

At variance with most kinetically constrained spin models, the spiral model is characterized by an equilibrium 
glass-jamming transition at a finite value of the control parameter: for p > p c an infinite cluster of frozen particles 
exists. In view of the more general issue of glassy systems, therefore, it is of a certain interest to study and understand 
the effect and the relevance of such a transition on the non-equilibrium dynamical properties. In so doing we uncover 
a quite interesting and rich scenario with unexpected far from equilibrium features. Actually, while the equilibrium 
dynamics (and the near-to-equilibrium one of the annealing process considered in Sec. IV Ap feels the presence of 
blocked regions the size of which diverges at the transition, leading to diverging relaxation times in the persistence 
function or in the autocorrelation, the non-equilibrium kinetics following an abrupt change of p, from below to above 
Pc takes a completely different route avoiding the percolative blocked configurations relevant in equilibrium. At the 
basis of this behavior is the reversibility property of any kinetics obeying detailed balance. Due to this fact, a blocked 
configuration cannot be connected dynamically with an unblocked one. As a consequence, a system prepared initially 
in any state without a spanning frozen cluster keeps evolving on a restricted ensemble of configurations without 
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such clusters. Evolution on a sub-ensemble of configurations (those with vanishing magnetization) is a key feature 
of coarsening systems. Actually, we found that the analogy between the spiral model and quenched ferromagnets 
goes beyond this. The non-equilibrium kinetics of the spiral model is actually a coarsening phenomenon where quasi- 
unidimcnsional strings of vacancies deform, diffuse and merge increasing their typical length L(t) in a power-law 
fashion (logarithmically in the extreme case p — ► 1). Observable quantities obey a dynamical scaling symmetry with 
L(t) acting as a scale factor, much in the same way as in usual coarsening. Interestingly enough, the non-equilibrium 
fluctuation dissipation plot is a non trivial (p-dependent) curve, qualitatively similar to the one observed in the Id 
Ising model. This feature may suggest that a similar mechanism is at work in the kinetics of the spiral model and of 
the Ising chain. It must be stressed, however, that in the former a non-trivial fluctuation-dissipation plot is observed 
in the whole glassy phase with p > p c , and not only in a limiting sector of the phase diagram (i.e. T — > 0) as in the 
Id Ising model, a feature resembling what is found in the Sherrington-Kirkpatrick spin-glass model [30j . 
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